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Abstract 

A model about the dynamic of vesicle membranes in incompressible viscous fluids 
is introduced. The system consists of the Navier-Stokes equations with an extra stress 
depending on the membrane, coupled with a Cahn-Hilliard phase-field equation in 3D 
domains. This problem has a time dissipative free energy which leads, in particular, 
to the existence of global in time weak solutions. The asymptotic behavior of these 
solutions is analyzed. Firstly, we prove that the w-limit set are critical points of the 
energy. Afterwards, by using a modified Lojasiewic-Simon’s inequality, we demonstrate 
the convergence of the whole trajectory to a single equilibrium. Finally, the asymptotic 
convergence of the phase field is improved. It is remarkable that all these arguments may 
be done without using strong regularity for the velocity and pressure variables. 

Keywords: Vesicle membranes, Navier-Stokes equations, Cahn-Hilliard equation, energy 
dissipation, convergence to equilibrium, Lojasiewicz-Simon’s inequalities. 

1 Introduction 

In this paper, we consider a model for the dynamic of vesicle membranes in incompressible 
viscous fluids, developed from the geometry of the membranes. This type of models was 
introduced by Helfrich [7]. The system consists of the Navier-Stokes equations with an extra 
stress depending on the membrane, coupled with a Cahn-Hilliard phase-field equation. 

‘This work has been partially financed by MINECO grant MTM2012-32325 (Spain) with the participation 
of FEDER 
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The membranes are formed by lipid bilayers. Under appropriate conditions, they with¬ 
draw into itself forming a sort of bag, named vesicle. The equilibrium configurations of vesicle 
membranes can be obtained minimizing the Helfrich bending elastic energy, under fixed sur¬ 
face area and volume constraints. 

Numerous studies have been devoted to this type of models and a detailed description of 
they can be seen in [1] and references therein. A phase function is also used to model the 
vesicle membrane as a diffuse interface in subsequent papers. In [9] and m, a coupled Allen- 
Cahn and Navier-Stokes problem is studied approaching both constrains, area and volume, in 
a penalized manner. In [T], a Cahn-Hilliard phase-field model is introduced, without taking 
into account the vesicle-fluid interaction. Now a Cahn-Hilliard-Navier-Stokes model will be 
considered. Since the volume constraint is implicitly satisfied in the Cahn-Hilliard equation, 
now only the surface area constraint must be approximated via penalization. 

We will prove that the resulting problem is thermodynamically consistent, because there 
exists a free energy (kinetic plus bending plus penalized one) which is dissipative in time along 
the trajectories. This fact is used to prove the existence of global in time weak solutions. 

On the other hand, the asymptotic behavior of the solutions is analyzed following the 
way of 0, m, m- We prove that the cu-limit set for weak solutions are critical points 
of the dissipative energy. After that, by using a modified Lojasiewic-Simon’s inequality we 
demonstrate the convergence of the whole trajectory to a single equilibrium. We consider 
of remarkable interest in this paper two facts: The study of the new model Navier-Stokes- 
Cahan-Hilliard modeling vesicles and the identification of a unique equilibrium point in the 
context of weak solutions. 

The current paper is organized as follows. We explain the model in Section 2 and give 
some preliminary results in Section 3. In Section 4, an energy equality and some global 
weak estimates are obtained that let us to prove the existence of weak solutions. Section 5 
is devoted to the study of convergence at infinite time for global weak solutions. We prove 
that the w-limit set are critical points. After that, by using a modified Lojasiewic-Simon’s 
inequality, we demonstrate that each trajectory converges to a single equilibrium. In section 
6, some global in time strong estimates for the phase are obtained, which let us to improve 
the convergence of the phase trajectory in a higher order space. 

2 The model 

We will analyze the case where the bending energy Ej, is given by a simplified elastic 
Willmore energy plus a penalization of the surface area constraint [5]: 
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( 1 ) 


Ebi(i>) = ^ [ dx + ]-M{A{4>) - of 

Jn £ ^ 

where F'{(p) = (|(/>P — !)</> for each (/> € H, being F{<:f>) = —— 1)^ the Ginzburg-Landau 
potential, M > 0 is a penalization constant, e is related to the interface width, and 

^(</>) = ^(||V0p + iF'(</.)) dx, 

approaches the surface area. 

Remark 1 The same results of this paper may be obtained replacing the surface area AiS) 
only by the first term A{fi) = [ ^\Vcj)\^ dx as in 0/. 

Jn 2 

By introducing the chemical potential 

SEkifi) 

^ ■ Sep 

we will introduce the following Navier-Stokes-Cahn-Hilliard equations in fl x (0,+oo): 


dtu+ {u-V)u— uAu— \wV(p -|- Vg = 0, 

( 2 ) 

V • M = 0, 

(3) 

dt(p + u ■ V(p — jAw = 0. 

(4) 


The coefficients > 0, A > 0 and 7 > 0 depend on viscosity, elasticity and mobility, respec¬ 
tively. The system ©-(HI) is completed with the boundary conditions 

'^\dU ~ 0) 9n(p\dQ = 0, dn^(p\dQ = 0, dnW\QQ = 0, (5) 

and the initial conditions 

u{0) = Uq, (/>( 0 ) = (pQ in n. ( 6 ) 

For compatibility, we will assume ito|an = 0 with V • rio = 0 and dn(pg\dn = 0. 

By integrating the rc-equation using the free-divergence S/ ■ u = 0, the non-slip 
condition u\qq = 0 , and the last boundary condition dnw\Qn = 0 , it is easy to deduce 

d f 

— / (p(x, t) dx = 0, 
dtJn ^ ^ 

hence one has the conservative character of (p in n, because the total volume is conserved: 

V{t) := / (p{x,t) dx= (po{x) dx := Vq € R. 

Jq J q. 
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On the one hand, for all ^ 

/6A{4>) - 


= / £Vcj)-V^ + -F'{cj))^, 


hence integrating by parts, if </> € and dncplon = 0 , one can identify 

d<p £ 

Note that, by using the boundary conditions for cj) given in ([5]), it also holds V/r • n|aQ = 0. 
On the other hand, for all 4>,(l> & H^, 
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M(-eA0 + -F'(0)^) + M{A{cl)) - a) i / eV(/. ■ , 


1 


hence, after some integrations by parts, using V(j)-n\QQ = 0 , V/r-n|gf 2 = 0 and Vcp-nlQfi = 0 , 
one can identify 

w := = - A/i + F'{(j)) + M {A{4>) - a)fj, 

o(p 


Remark 2 The variational derivatives w := 
L'^(Q.)-functions via the L^( 0 ) scalar product. 

In particular, we can decompose 


6(f 


and p, := 


5A{4>) 


have been identified as 


w = £/S?(j) + G{(fi) 

with 


(7) 


G{fi) := --/^F’ifi) + \piF'{fi) + M{A{fi) - a)pi. ( 8 ) 

e £^ 

Since dnfifin = 0, in particular = 0, hence 

[ -AFffi) dx= [ -F"{fi){Vfi ■n)dx = 0. 

Jo. Jan 

Therefore, integrating ([7|) and ([8]), 

I w dx = I G{4>) dx = — f p. F'(fi) dx + M{A{(j)) — a) f fidx. (9) 

Jn Jn ^ Jn Jn 

By using w = eA^</> + G{(f) as auxiliary variable, we can rewrite the problem ([2])-([51) as 

dtu+ {u-S/)u — oAu — XivS/fi + Vq = 0, (10) 

V-M = 0, (11) 

dtfi + u ■ \/(f> — jAw = 0, ( 12 ) 

£A‘^(j) + G((/)) — rc = 0, (13) 
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With the aim to use the conservation of the problem, if we denote mo = := tvw / 4>o{x) dx, 

l“l Jn 

we introduce the following mean-value variables for the phase-field problem: 


^{x,t) := (f>(x,t) — rriQ and z := w — {G{(p)). 

From ([9]), {G{4>)) = (fx F'((f>)) +M{A((I)) — a){ij). Reciprocally, given ('ip,z) we can recover 
{ip, w) as V’ = V’ + "^0 and w = z + {G{(j))). 

By rewriting the equations (I10l) - (ll3p in these new variables {ip, z) we arrive at the problem 

dtu+{u-V)u —v/S.u—\zVip+ Vq = (), (14) 

V-ii = 0, (15) 

dtip + u-Vip — ^/S.z = ^, (16) 

e/SP^ip -|- G{ip) — z = 0, (17) 

where 

q = q + X{G{(p))ip 

and 

G{ip) := G{ip + mo) - {G{ip + mo)), 
completed with the boundary and initial conditions 


”*^190 — d) dn'^\dQ. — 0, dnAi.1p\QQ — 0, dn^\dQ, — 0) 

m( 0) = ito, ■0(0) = ipo := 00 - (0o) in 11. 


Observe that, J ip dx = 0 and / z dx = 0. 
Finally, by (denoting 


/ 

Jn 


then, 


z = 


Eh{ip) = Eb{ip + mo) 

6Eb{iP) 


dip 


= eA^iP + G{iP), 


where the identification (j21ll has been done via the L^(ll)-scalar product. 


(18) 

(19) 


( 20 ) 

( 21 ) 


3 Some preliminary results 

Let us to introduce the following notations; 

• In general, the notation will be abridged. We set = L^’(ll), p > 1, = E[q{Q), etc. 

If X = X(ll) is a space of functions defined in the open set 0, we denote by LP{0,T; X) 
the Banach space LP{0,T; X{f})). Also, boldface letters will be used for vectorial spaces, 
for instance = L‘^{Q)^. 
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• The L^-norm is denoted by | • |p, 1 < p < oo, the i7™'-norm by || • \\m (in particular 
I • I2 = II • llo)- The inner product of L‘^{Q,) is denoted by (•, •). The boundary H^{dQ)- 
norm is denoted by || • ||s;ao. 

• We set V the space formed by all fields u G satisfying V • u = 0. We denote H 

(respectively V) the closure of V in (respectively H^). H and V are Hilbert spaces 
for the norms | • I 2 and || • ||i, respectively. Furthermore, 

i? = {m G V-M = 0, M-n = 0on 9^}, {ti G H^; V • u = 0, m = 0 on 9^} 


• From now on, (7 > 0 will denote different constants, depending only on the fixed data 
of the problem. 

Let us consider the following mean-value spaces: 


Ll = e j u; = o|, 

e J w = 0^ k>l, 


Hf = {w e i7*(ri); dnw\dn = O} 

= {w e iL*; dnw\dn = 0, dnAw\dn = 0} k = 3,4, 5,6. 

We will assume H regular enough to use the regularity of the two following elliptic Laplacian- 
Neuman and Bilaplacian-Dirichlet-Neumann problems, respectively: 

A^V’ = / ill 


—Az = f in n 
dnZ\dQ = 0, [ z dx = 0, 


— 0 ) dnAlplQQ — 0 , 


/ 

*. 1/ Q 


Ip dx = 0 


where / G L^(n). From the Ff^-regularity of the first problem, the following norms are 
equivalents: 

||z ||2 ~ |A 2;|2 in Hf, (22) 

and from the H^, and Ff^-regularity of the second problem, 

||V^||4 ~ |AV |2 iniLl, IIV^IIs Ri IIA^V^IIi iniLf, HV’He ~ ||A 2 V :||2 in iL|. 

Now we study the elliptic problem relating fj, to V’: 

Lemma 3 Given fi, we consider V’ as the solution of problem 

—eAif) + -|- mo) = ^ in ^ 


(23) 


dn'iPldn = 0, / dx = 0. 
Jn 
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Then, 


\\ip\\i < C{1 + 1^1 

2 < <^(1 + \ t \2 + IIV’lll); 


(24) 

(25) 


Proof. Firstly, by taking -0 + uiq as test function, we obtain 

e\V'ip\l + + + mo) + + 

By using Young, Holder inequalities and the Poincare inequality for mean-value functions we 
obtain: 

C'^IIV’II? + + molt < + (^ + ^)IV’ + moll < + ^IV’ + "iol4 + C, 

hence, 

IIV'lli < c\fj,\l + c 

and (1241) holds. Secondly, from the regularity of the problem (bootstrap’s argument) 

1 

—sA'ijj = -+ mo) in 12 

= 0 > / dx = 0, 

Jn 

we obtain that ||V ’||2 < C'd/^b + + mo)\ 2 )- From the dehnition of F', 

|F'(V’ + mo)|2 < C'dV’ + "lole + IV’ + "lob) < C'(||V’ + "io||? + IV’ + "I 0 I 2 ) < <^(1 + 

Therefore, ||V’||2 < C{^ + |m |2 + HV’lli) arid (f25]l holds. 


We introduce a modihed Lojasiewicz-Simon’s inequality. 

Lemma 4 (Lojasiewicz-Simon inequality) LetS be the following set of equilibrium points 
related to the bending energy given in [W\) : 

<S = {V’ G H 2 {TI) : eA^V’ + G'(V’) = 0 a.e in 12}. 

Let V’ G 5 and K > 0 fixed. Then, for any two sufficiently small constants (3 and 5, there 
exists C > 0 and 9 G (0,1/2) (depending on ip, ft and 5), such that for all ip G Ffl 
llV’lb < K, \\ip — ip\\i < fi and \Ei,{ip) — Ei,{ip)\ < 6, it holds 

\Ef,iiP)-E,ilp)\^-^ <C\z \2 (26) 

where z = z{ip) := eAfip + G{ip). 
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Proof. 


Step 1 : For /?! small enough, there exists C > 0 and 9 € (0,1/2) (depending on ip and (3i) 
such that for all p! G with \\ip — V'lU < /^i? then (f^) holds. 

The proof of this step is based in the generic result Theorem 4.2 of [ 8 ] for the following 
spaces and operators (with the same notation that in [ 8 ]): 

H = X = Ll{n), X = Ha = Hi, 

A = A'^ e X AiP = AV g H and (V’, Oa = A'^^)l 2 V’ € D{A), 

G A i-G S{ip) = Ei){'ip) G M, such that 

Ebii^) = Ebipj + mo) = ^ [ (-eAV’ + -T'(V’ + mo))^ dx + ^-MiAipj + mo) - a)^ 

2 e Jq e 2 

M=S' -.ipeXe^ M{ip) = eA^ + 7^(V’) - {n{pj)) G H, 
where TZ{'ijj) = G{'ip + mo) hence G(V’) = “ (^('0))- Finally, = M.'i'ip), where 

M'ip ,): G A ^ M'im) ■= eA^e+ n'im) - (H'im)) e h, 

with 

H'WiO = -^[F'”{^P + mo)VV’ • Ve + 3|VV^|2^ + F"{xP + mo)A^ + F'"(V’ + mo)(AV^)e] 

+ ^[(F"(V’ + mo))^ + F"'(V’ + mo)F'{'ip + mo)]^ 

If 1 

+M(—eA'i/’H— F'{ip + mo)) / VV'• dx + M(A(?/i + mo) — q;)(— eA^ H—F"(V’ + mo)0- 
^ Jn ^ 

Note that Xi^tp) is a Fredholm operator of index zero, because Ai^ip) is the sum of the 
invertible operator A and a compact operator. 

Moreover, the map TZ : ip ^ X t-A- M'{ip)A~^ G C{H) is well-posed because A~^ G C{H-,X) 
and Xi'ptp) G C{X;H). It remains to be proved that TZ is (sequentially) continuous. Let 
—>■ 'i/ in A as n —>■ oo. Then, 

\\n{P;^)-nmc(H) = \\M'{P^n)A-^ - M’{iP)A-^ciH) 

< \\M'{P;n) - M'm\ciX;H)\\A-^\\ciH;X) 


8 


Since 


7e'(V’n)(0 - + mo)V(t/>„ - V^) + (F'"(V>n + mo) - F'"(V’ + mo))VV^]Ve 

-^[3(VV^„ + V^P){V^Pn - VV^)e + (F"(V'n + mo) - F"(V^ + mo))Ae] 

-+ mo)A('!/’n - V’) + (V’n “ V’)A'!/’]? + ^[-F"(V' + mo)3((?/)„ + mo)^ - (V’ + mo)^) 

+(F"(V>n + mo) - F"{'ip + mo))(3(V'n + mo) - 1) + Q{i^n + mo){F'{ipn + mo) - -F'(V’ + mo)) 

+6{'^|Jn-'lp)F'{^l; + mo)]^ + M[{-£A^lJn + -F\‘^|; + mo)) [ SI{ipn - i dx 

£ Jn 

{-eAiif^n -'!/’) + -{F'itpn + mo) - F'i'tjj + mo))) [ SI'ipSI^ dx] 

e Jq 

+M{-eA^+-F"iiP + mo)0 [ - VV’)(VV’n + VV’) dx 

£ Jn ^ 

+M{-sA^ + -F'^i/j + mo)0 / i-((V'n + mo)^ + (V' + mo)^ - 2)(V'„ + ^/; + 2mo)('0n, - 'ip) dx 

e JO 4e 

+-(A(V' + mo) - a){F"{'ipn + mo) - + mo)) 


and the differences F'{ipn + mo) — F'(■;/’ + mo), F''{ipn + mo) — + mo) and F"'{'ipn + 

mo) — F'"{'ip + mo) contain the factor ipn — 'P’, then 


WM'i'ipn) - M'{'ip)\\c(x-,H) = sup 

Cex\{ 0 } 


||Af'(V^.)(0->('(V’)(0k 


lien 


X 


= sup 
$GX\{ 0 } 


|g(^n)(e)-gw(e)i2 
iieik 


< C'dIV'n - V'IIlZ + II^An - V'IIhI + llV'n “ 'IPW^)- 


Therefore, ||A^'(V’n) — ■M'{'4’)\\c{x-,h) — >■ 0 as n ^ oo if V’n — t V' in (even ii ipn ^ 'ip in 
H^), hence the continuity of the operator F is deduced. 

For any ip € FK^i) satisfying ||'0 — 'i /)||4 < /3i, owing to Theorem 4.2 of [8], there exists 
C > 0 and 6 G (0,1/2) such that: 


\E,i^p) - <c\\£'m\H = c\zm 2 


and (IU2]) holds. 

Step 2: (Relaxing the local approximation \\ip — ^p\\i < /3i by llV'lb < K, Wp; — ^p\\l < pi and 
\EbiP’) — Eb{^p)\ < 6). There exits C > 0 and 9 € (0,1/2) (depending on ip, pi, K and 6) such 
that if ip € with HV'Hs < K, Wp; — ?/||i < jd and \Eb{'p) — Eb{'ip)\ < 5, then (f62)) holds. 

In this step, we follows Lemma 4.4 of m but imposing the “proximity” condition between 
Ip and ^p only in the iL^-norm instead of the F^-norm as in m by using the convergence of 
the energies. 

Firstly, there is a constant M > 0, such that 

IIV^ - 'ipWi < M|A^(V^ - '!/;)|2 
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(here, J ijj dx = J “ip dx has been used), hence, 

\z{ip)\2 > \eA‘^{ip -^p)\- \G{ip) - G{ip)\2 >^\\'ip- 'iplU - IG(V’) - G(V’)|2. (27) 

Secondly, we are going to bound \G{ijj) — G{ijj)\ 2 - For this, 

G(V’) - G{lp) = -klF’^pj + mo)V(V' + ^)VtP + 6V'(V^)2] 

—2k[F"{'i/j + mo) Alp + 3{ip + pj + 2mo)'ipA'ip] 

k — — 

+^[F"{ip + mo)ip{{ip + mo)^ + (V’ + mo)i'ip + mo) + (ip + mo)^ - 1) 

+3ip{ip + ip + 2mo)F'{ip + mo)] 

Me^ 1 , C — 

H— —{-A'lpF-F^'ip + mo)) / V{-ip +'ip)Vip dx 

2 £ Jq 

+M£{B{'ip) - /3){-Aip + ^{F'{p; + mo) - F^p; + mo))). 

Then, there is a constant Gi = Gi{K) > 0 such that 


|G(V’)-G(^)|2<Ci(i7)||V>-V’l|2. 


In particular, interpolating Hf between Hi and 77|, 

\G{p;) - G{p;)\2 < Ci{K) \\p; - ~ ^Gi{K)‘^\\p; - p;\\i. 


Let f3i > 0, 9 € (0,1/2) given in Step 1, by choosing (5 > 0 and P > 0, both sufficiently 
small, such that 




and 


rl-0 


< 


£/3i 

4M’ 


then, for any pj G H^i^kl) satisfying \\pj — pj \\2 < P and \Ei,{pj) — Eh{p>)\ < 6 we have 


\GiP;)-GiP;)\2< 


2M' 


- V'lU + 


ePi 

AM 


(28) 


and 

\E,{P;) -E,ilp)\^-<^ < ^. (29) 

There are only two possibilities: either ||'!/; — '0||4 < Pi and then (l62|) holds by using Step 1; 
or Wp; — '0114 > Pi- In this latter case, from ([?7|1 and (1^ 

k(0)|2 ^ ^IIV’ - V’lk - ^C'l(iL)2||0 - 0||l 
^ 2M AM ~ AM' 


On the other hand, from (1^ . > \Eh{pj) — Eh{p>)\^ ^ = \£{pj) — £"(0)1^ ^ hence, (f6^ 

holds. ■ 
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4 Weak Solutions 


Definition 5 Let uq € H and ijjo = 4^0 ~ S Hf, we say that {u,'ip,z) is a global weak 


solution of (fT^ - (fTn|) in (0,+cx)), if 

ue L‘^{0,+oo;V)riL°°{0,+oo-,II), z e L‘^{0,+ 00 ] Hi) 

G L°^{0,+oo;Hf), 

satisfying the variational formulation 

{dtU,u) + {{u-\/)u,u) + i'{Vu,Vu) — X{z\/'if,u) = 0 y u ^ V, (31) 

{dtip,z) + (u-Vif,^) + j{Vz,Vz) = 0, yzeHl (32) 

siAif, Af)) + (G(V'), V') -{z,'if) = 0, y-ip e Hf, (33) 

and the initial conditions m- 


Observe that, from (1301) . (I3TI) and ([32]) . one has dtu ^ L^/^(0, + 00 ; V') and dtfj G L‘^{0, + 00 ; (Hi)'), 
hence in particular the initial conditions ()19[) have sense. 


4.1 Energy equality and large-time estimates 

In a formal manner, we assume that (tt, ip, z) is a regular enough solution of (|14p - (|19p . By 
taking u= u, z = z and fp = dtip as test function in (|3Tp . (f32p and (l33p respectively, 

- X{zV'ip,u) = 0, 

{dtip, z) + {u- VV’, z) + 7I Vz|2 = 0, 

+ {G{ip),dt'ip) - {z,dt^p) = 0. 

Adding the hrst equality plus the second and third ones multiplied by A, the term (z, dt^p) 
cancels, as well as the nonlinear convective term {u - Vip, z) with the elastic term —(z VV’, u), 
arriving at the following equality 

+ V|A,/.g) + X(a()p),8,i>) + v\Vu(t)\l + A7|V2(i)|| = 0. (34) 


Si„ce 

0(p 


z, in particular 


^Ebiipit)) = {^-^^^dt'ip) = (z^dtip) = s~\A4;\l + {G{ip),dt'ip). 
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We define the total free energy E{u, ijj) = Ek{u) + XEb{xjj), where is the bending energy 

defined in ([20]) and Ek{u) = 5 / \u\‘^ the kinetic energy. Then, equality ([M|l can be rewriten 

Jn 

as the following energy equality: 

^E{u{t),'ijj{t)) + n\Vu{t)\l + X-f\Vz{t)\l = 0, (35) 

which shows the dissipative character of the model with respect to the total free energy 
E{u,'ip). Moreover, assuming the initial estimates (tto,'0o) in i?x Hi, the following uniform 
“weak” bounds in the infinite time interval ( 0 , + 00 ) hold: 

uin L°°{0,+oo-,H)nL^{0,+oo-,V), in L^{0,+oo-, Hf), z in L^{0,+oo-,H^). (36) 

4.2 Additional estimates for ip in iff 

The idea is to use previous weak estimates ([36l) in the '0-equation m- Since 0 G 
L°°(0,+ 00 ; 4t^), in particular 0 G L°°(0,+ 00 ; L°°), then F'{ip), E'{ip) and E"'{ip) are also 
bounded in L°°(0, + 00 ; L°°). Therefore, we have 


|G'(0)|2<C, 

(37) 

|VG(0)|2<C(1 + ||0||3), 

(38) 

|AG(0)|2<C(1 + ||0||4). 

(39) 


From the 0-equation (fTTjl . by using ([23|), (f37|l . (I38|) and (f39]l . we obtain 

II0II4 < C {1 + |G(0)|2 + \z\2) < c(l + |2:|2)- 

In particular, from (|39ll and Poincare inequality for the mean-value function z, 

|AG(0)|2 < C (1 + |z| 2 ) < C{1 + |Vz| 2 ) (40) 

On the other hand, again from (fTTll . by using (l38l) and the interpolation inequality ||0||3 < 
||0||2^^||0|iy^, we deduce ||0||5 < 0(112:111-|-||0||2'^^ 11011^^) and, therefore, since ||0||2 < C and 
II 0 II 4 < 0(1 -|- | 2 ;| 2 ), then ||0||5 < 0(||2;||i -|- 1). In particular, 

0GLL(O,+oo;F5). (41) 

For instance, weak solutions furnished by a limit of Galerkin approximate solutions satisfy 
the corresponding energy inequality related to (155]) (changing in ([551) the equality = 0 by < 0) 
and this inequality energy suffices to prove rigorously all previous estimates ([36]) and ([4X1) for 
the Galerking approximations. 

Gonsequently, fixed the initial datum (tiO)0o) G i? x Hf, by using a Galerkin Method 
and proceeding in analogous way to Subsection 3.3 of [2], one can prove existence of weak 
solutions of dH-dlS]) in ( 0 , -|-oo). 
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5 Convergence at infinite time. 

From the energy inequality (1351) . we have 

E{u{t),'ilj{t)) - E{uo,'iIjo) + f {i/\Vu\l + X^\Vz\ 2 ) dr <0, Vt>0. (42) 

Jo 

Therefore, there exists a number E^^ > 0 such that the total energy satisfies 

E{u{t),il;{t)) \ Eoo R astj+oo. (43) 

The w-limit set of a fixed global weak solution, {u, ijj), associated to the initial data (tfO) V’o) £ 
Hx Hf, is defined as follows: 

Uj{u,xjj) = {{Uoo,'4’oo) € Hx Hf : t Too s.t. 

{u{tn),'ip{tn)) {Uocipoo) Weakly in H x Hi). 

Let S be the set of critical points of the energy E{u,'tp): 

S = {(0, Ip) : Ip ^ (^) • + G{'ip) = 0 a.e in 0}. 


Theorem 6 Assume that {uo,'ipo) £ Hx Hf. Fixed {u{t),ip{t)^ z{t)) a weak solution of m- 
m in (0,+oo), then, uj{u,ip) is nonempty and Lo{u,ip) C S. Moreover, for any {Q,ip) € S 
such that (Oj'ip) G a;(rt, V'), then Eb{'ip) = E^o holds. (In particular, u{t) —>■ 0 weakly in T^(fl) 
and Eb{'ip{t)) — Eb{'p) in IR as 11 +ooj- 

Proof. Since E{u{tn),'ip{tn)) < E{uq,'iPq), it exists at least a (non relabeled) subsequence of 
tn, such that 

u{tn) —> Uoo weakly in H, ipptn) —>■ ipoo weakly in (44) 

for suitable limit functions (wooV’oo) £ Hx Hf. Let us consider the initial and boundary 
value problem associated to (flDl - lfTOll on the time interval [tn,tn + 1] with initial values u{tn) 
and 'tp{tn). Setting 


Un{t) := u{t + tn) -ipnit) ■■= ip{t + tn) and Zn{t) := ziipit + tn)) for a.e. t G [0,1], 

then, {un,'ipn) is a weak solution to the problem (fT4)l - (fT^ on the time interval [0,1]. In 
particular, from (|42]l . 




/ {v\V Un\l + \-i\VZn\l) dr 

Jo 

{v\Vu\^ + AqlVzjl) dr < E{tn + 1) — E{tn) —>-0 as n —>■ oo. 


(45) 
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In particular, Vua{t) —>■ 0 strongly in L^(0, and Vz(tn) —>■ 0 strongly in L^(0, 

whence, by Poincare inequality, 

—>■ 0 strongly in L^(0,1 ; V) and Zn(t) —?> 0 strongly in L^(0,1; HI). (46) 

On the other hand, from the large time estimates (I36p . ipn is bounded in L°°(0,1; Hf) and 
is bounded in -L^(0,1; H^)- Then, from the V'-equation (IIZl), dt'ipn is bounded in L^(0,1; (Hi)'), 
therefore, by using the Aubins-Lions compactness theorem, 

ipnit) Ip strongly in C°([0,1]; i?i). 

In particular, '0n(O) —)• V^(0) in H^, hence 'tp{0) = ipoo is attained because V’n(O) = ip{tn) —^ V’oo 
in Pf^-weak. Moreover, dtipn converges weakly to dti) = 0 in T^(0,1; (Hi)'), therefore, 'ip is 
constant in [0,1]. Consequently, 

= ^poo for all t G [0,1]. (47) 

Finally, by using convergences (I4U1) and G{'ipn) —>■ G{ip) weakly in T^(0,1;T^), taking limit as 
n —)■ 0 in 

{Zn,Pp) = {Alpri, A'^) + (G(V’n), ^), G 

we deduce that (A'0„, A'0) converges to (A'0, A^/i) and also, 

{Alp, App) + {G{'ipn),Pp) = 0, G Hf, a.e. t G (0,1), 

that is {A'lpoo, Ap)) + {G{ipoo),'ip) = 0, 'i^p G Hf. Therefore, ipoo is a weak solution of the 
elliptic problem 

( A^p) + G(V’) =0 in n 

I SnV’jan = 0, dnAp;\QQ = 0, j pj dx = 0. 

Finally, by using a bootstrap’s argument, the regularity p^oo G is deduced and the proof 
is hnished. ■ 

Theorem 7 Under the hypotheses of Theorem 0 there exists a unique limit pj G such 
that p}{t) pj in weakly as t p +oo, i.e. io{u, p)) = {(0, pj)}. 

Proof. Let (0, p;) G u}{u, p;) C S, i.e, there exists tn p +oo such that uptn) —>■ 0 weakly in 
and p’{tn) —>■ V’ fo Hf weakly. 

Without loss of generality, it can be assumed that E{u{t),pj{t)) > E{0,pj){= E^o) for all 
t > 0, because otherwise, if it exists some t > 0 such that E{u{t),p^{t)) = E{0,pj), then the 
energy inequality (t35|) implies 

E{u{t),p’{t)) = E{0,p;), |V'u(t )|2 = 0 and |Vz(t )|2 = 0, for each t > t. 
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Therefore, u{t) = 0 and z{t) is constant Vt > t. In particular, by using the z-equation, 
dtipit) = 0, and hence V’(^) = 4’ for each t > t. In this situation the convergence of the 
V^-trajectory is trivial. 

Assuming E(u(t),4j{t)) > E{0,'ip){= E^o) for all t > 0, the proof is divided into three 
steps. 

Step 1: Let 4j € Hf such that \\4>{t) — ^p\\l < j3, for each t > ti > 0, where (3 > 0 is the 
constant appearing in LemmaO (of Lojasiewicz-Simon’s type), then the following inequalities 
hold: 

j^(^(E{u(t),m) -mW) +ce (|Vt/(t )|2 + |Vz(t)| 2 ) <0, Vt G (ti,oo) (48) 

\\dt'i)\\iHlY<j(E{u{ti),ip{ti))-E{0,il;))f, Vt2>ii, (49) 

where 9 G (0,1/2] is the constant appearing in Lemma [H 
Indeed, the energy inequality (f35]l can be written as 

j^{E{u{t),m) -Eoo) + C {\Vu{t)\l + \S/z(t)\l) < 0, 
hence, in particular, from Poincare inequality 

^{E{u{t),'ip(t)) - Eoo) +C{\u{t )\2 + |^(t)| 2 ) {\Vu{t )\2 + |V 2 ;(t)| 2 ) <0, Vt > 0. 
Therefore, by taking the time derivative of the (strictly positive) function 

H{t) := (E{u{t),4;{t)) - E^f > 0, 


we obtain 

^^+9(E{u{t),m) - E^f-^C{\u{t )\2 + \z{t)\ 2 ) (|Vn(t )|2 + \Vz{t)\ 2 ) < 0, Vt > 0. (50) 


On the other hand, since the unique critical point of the kinetic energy is m = 0, by taking 
into account that \Ek{u) — L'fc(0)| = that there exists a constant K such that 

|w(t )|2 < K, then 

\Ek{u{t)) - .Lfc(0)|^“'^ = = ^|u(t)l 2 “^V(i)| < C\u{t )\2 Vt > 0. 


This estimate together the Lojasiewicz-Simon inequality (given in Lemma d]): 


\Eb{4^{t)) - Eb{4j)\^ ^ <C{\u{t )\2 + \z{t)\ 2 ), Vt > ti. 
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give 


{E{u{t),m) - Eoof-^ > \Ek{u{t)) - Ek{0)f-^ + \Ebim) - EbWf-^ 
>C{\u{t)\2 + \z{t)\2)~'^ 


(51) 


Applying (fSTj) in (f50]) . 


dH{t) 

dt 


+ ec{\vu{t)\2 + \vzit)\2) <0, vt>ti 


and (H5D is proved. 

Integrating (H8l) into for any t 2 > h, we have 

rt2 


{E{u{t 2 )Mt 2 ))-E^f + ec l\\Vu{t )\2 + \Vz{t)\ 2 )dt<{E{u{ti),'il^{ti))-E^f. (52) 

Jti 

From the z-equation, dtip = —tt-VV’ +Az, and by using the weak estimate ||i/'(t )||2 < C, then 


\m\\^Hiy<C{\Vu\2 + \Vz\2) 

By using this inequality in ([5^ . then is attained. 

Step 2 : There exists a sufficiently large no such that ||V'(t) — '0||i < (3 and |n(t )|2 < K for 
all t > trio (1^ E given in Lemma Up. 

The bound |n(t )|2 < K Vt > 0 has been already obtained from weak estimates of u in 
L°°(0, +oo; IT). We now focus on the bound for ||i/)(t) —Since 'ip{tn) —in (strongly) 
and E{u{tn),tp{tn)) —>■ E^o = Eh^tp), then for any e € (0,/3), there exists an integer N{e) 
such that, for all n > N{e), 

Wi’itn) - tpWi < e and ^{Eb{u{tn),'tp{tn)) - E^o)^ < e. (53) 

For each n > N{e), we define 

tn := sup{t :t>tn, IIV’(s) - V’lli <13 Vs E [tn, t)}. 


It suffices to prove that = +oo for some uq. Assume by contradiction that tn <tn < +oo 
for all n, hence ||V’(in) — '0||i = /3 and \\ip{t) — ip\\i < (3 for all t E \tn,tn)- By applying Step I 
for all t E [tn-,tn]-, from and (|53|) we obtain, 


Therefore, 



Wdi-ipW^Hiy < Ce, 


Mny N{e). 


_ _ rt„ 

U(tn) - V'll(Hl)' < U{tn) - V'll(Hl)' + / WdtipWiHiy < (1 + C)e, 

J tn 
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which implies that lim„^+oo ||'0(^n) — = 0- Since V’ is bounded in L°°(0,+oo; 

{'ip{t))t>t* is relatively compact in H^. Therefore, there exists a subsequence of il^(in), which 
is still denoted as ipitn), that converges to ijj in Hence, \\ip{tn) — V’lli < /? for a sufficiently 
large n, which contradicts the definition of tn- 

Step 3 : There exists a unique ij: such that 'ip{t) — >■ 'll; weakly in as +oo. 

By using Steps 1 and 2, can by applied, for all H > to > Ino, hence 

mti) - i^{to)\\(Hl)' < ! IIStV'llcHi)'^ 0, asfoTi^+oo- 

Jto 

Therefore, (V’(f))t>4no is a Cauchy sequence in {Hi)' as f t +oo, hence the (-ff* )'-convergence 
of the whole trajectory is deduced, i.e. there exists a unique V’ £ {Hi)' such that V'(t) —)• V’ i^i 
{Hi)' as f t + 00 . Finally, the weak Ff^-convergence by sequences of 'ip{t) proved in Theorem^ 
yields 'ip{t) ijj in H^ weakly. ■ 


6 Higher estimates only for the phase variable 

6.1 Global in time strong estimates for Tp 

By adding the z-equation (l3^ tested by dll’ G H\, and the V'-equation (j3^ by G Hf 

(see (fTSjl b integrating twice by parts the term {G{p) — z, Adtl’), taking into account that 
Vdtp ■ n\Qn = 0, VG{'p) ■ n|ao = 0 and Vz ■ n|ao = 0, then 

(G(V’) - 2, Adtip) = (AG(V^) - Az, dtip). 

Then, the term {Az, dtp) cancels, remaining: 

e^^lVAV^ll + \dtP\l = -{u-S/p, dtp) - {AG{p),dtp). 

In particular, 

e^lVAV’li + \dtp\l <G{\u- Vp\l + \AG{p)\l). (54) 

We bound the convective term as 

l«-VV’li<|M|^|VV’|i<C||n||2. (55) 

From (flT)) . ((23]) . (f39]l and (|55]) . we have that HV’He < G{1 + 11^^114 + \dtp \2 + ||“||i)- By 
using the interpolation 11^^114 < G||V^||2^^||'0||y^||V'||g^‘^, one has 11^^114 < G||V^||g^^. Therefore, 

llV’lle — ^(1 + + ll^lli)- (b6) 
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By using (1^ in (l5H) . 


A 

dt 


iVAV^li + CollV’ll^ < c (1 + ||n||2 + 1^. VV^li + \AGml) 


and owing to (P0]l and (1551) . we have 

d 


dt 


\VA^l;\l + Co\ml<C{l + \\u\\l + \Vz\l). 


( 57 ) 


Since ||V^||3 is equivalent to (|AV'|2 + |VAV'|2) and, taking into account the weak estimates of 
ijj, llV’lla is equivalent to (1 + |VAV^|2). Then, from ([57)1 


A 

dt 


3 + C'oll'fAlle < 17 (l + IImII^ + IV^I^) . 


(58) 


By denoting 


B(t):=||rx||? + |Vz|i, 


then (1581) is rewritten as 


$' + Co^> < C{l + B). 

Multiplying (I59p by and integrating in time, we obtain 


(59) 


$(t) < $(0)e-^°^ + f e^“*(l + B{s)) ds. 

Jo 

In particular, <I>(t) < <I>(0) + (7(1 — + C f B{s)ds. Since B{t) € L^(0,+oo), we have 

Jo 

that G L°°([0, Too)). Moreover, integrating (fSTjl and ([^ in [0,t], we obtain 

V^€L-(0,+oo;77f), € LL(0,+cx); F|), G LL(0, +oo; L^). (60) 

Finally, from ([T7)l . 

^ € 7/;^oc(0) Too; Hf). 

In particular, using this improved estimates for the phase, the phase equation is satisfied 
point-wisely a.e. t G (0,+oo), if the data V’o is sufficiently regular. 

Remark 8 For this model, it has been possible to obtain higher estimates for the phase with¬ 
out improving estimates for the velocity and pressure. 


6.2 Improving the convergence of the phase trajectory 

The previous extra regularity obtained for if allows to obtain a more regular equilibrium 
point by using a Lojasiewicz-Simon inequality more demanding than Lemma 01 where the 
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proximity hypothesis between and ip is imposed in Hf. Fixed the initial data {uQ,ipo) G 
H X H2, the cj-limit set and the set of equilibrium points are defined by: 

Uj{u, Ip) = {{Uoc,ipoo) e Hx f +00 s.t. 

(blj 

{u{tn),ip{tn)) iuoo,ipoo) Weakly in x 
S = {(0, Ip) : Ip ^ -f^2 (^) • ^kA^ip + G{ip) = 0 a.e in 17}. 


Lemma 9 (Lojasiewicz-Simon inequality) If (0, ip) G S, there are three positive con¬ 
stants C, P, and 6 G (0,1/2) depending on ip, such that for all ip G and \\ip — ip \\2 < j3, 
then 

\E,{iP)-E,{fp)\^-^ <C\z \2 (62) 

where z = z{ip) := £/S?ip + G{ip). 

This result is rather classical. In fact, a similar Lojasiewicz-Simon’s lemma is proved in 
Lemma 5.2 of [T^ . 

In this setting, theorems [6] and [7] can be extended in the following way: 

Theorem 10 The set uj{u,ip) given in i61\) is nonempty and uj{u,ip) C S. Moreover, for 
any {0,ip) G S such that {0,ip) G uj{u,ip), then Eh{ip) = E^^ holds. 

Theorem 11 There exists a unique ip such that {0,ip) G S and ip{t) ip in weakly as 
11 +00, i.e. uj{u, Ip) = {(0, Ip)}. 

7 Conclusions and Prospects 

For the Navier-Stokes-Cahn-Hilliard model introduced in this paper we have proved the 
convergence of the each trajectory to a single equilibrium point for global weak solutions. 
Moreover, the regularity for the phase is improved from the energy inequality without the 
need of more regularity for the velocity and pressure variables, and hence large time or large 
viscosity constraints are now unnecessary. 

Starting of the results obtained in this paper, it seems achievable to obtain (rational) 
convergence rate estimates of the convergence of trajectories in a similar way to [6]. Finally, 
it would be interesting to study if local minimizers of the elastic bending energy are the 
stables, as has been done in [12] for a Navier-Stokes-Allen-Cahn problem modeling vesicle 
membranes. 
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